36M-pixel synchrotron radiation micro-CT for whole secondary pulmonary lobule visualization from a large human lung specimen

A micro-CT system was developed using a 36M-pixel digital single-lens reflex camera as a cost-effective mode for large human lung specimen imaging. Scientific grade cameras used for biomedical x-ray imaging are much more expensive than consumer-grade cameras. During the past decade, advances in image sensor technology for consumer appliances have spurred the development of biomedical x-ray imaging systems using commercial digital single-lens reflex cameras fitted with high megapixel CMOS image sensors. This micro-CT system is highly specialized for visualizing whole secondary pulmonary lobules in a large human lung specimen. The secondary pulmonary lobule, a fundamental unit of the lung structure, reproduces the lung in miniature. The lung specimen is set in an acrylic cylindrical case of 36 mm diameter and 40 mm height. A field of view (FOV) of the micro-CT is 40.6 mm wide × 15.1 mm high with 3.07 μm pixel size using offset CT scanning for enlargement of the FOV. We constructed a 13,220 × 13,220 × 4912 voxel image with 3.07 μm isotropic voxel size for three-dimensional visualization of the whole secondary pulmonary lobule. Furthermore, synchrotron radiation has proved to be a powerful high-resolution imaging tool. This micro-CT system using a single-lens reflex camera and synchrotron radiation provides practical benefits of high-resolution and wide-field performance, but at low cost

Anonymization server system for DICOM images

We have developed an anonymization system for DICOM images. It requires consent from the patient to use the DICOM images for research or education. However, providing the DICOM image to the other facilities is not safe because it contains a lot of personal data. Our system is a server that provides anonymization service of DICOM images for users in the facility. The distinctive features of the system are, input interface, flexible anonymization policy, and automatic body part identification. In the first feature, we can use the anonymization service on the existing DICOM workstations. In the second feature, we can select a best policy fitting for the Protection of personal data that is ruled by each medical facility. In the third feature, we can identify the body parts that are included in the input image set, even if the set lacks the body part tag in DICOM header. We installed the system for the first time to a hospital in December 2005. Currently, the system is working in other four facilities. In this paper we describe the system and how it works

Computer aided diagnosis for severity assessment of pneumoconiosis using CT images

240,000 participants have a screening for diagnosis of pneumoconiosis every year in Japan. Radiograph is used for staging of severity in pneumoconiosis worldwide. This paper presents a method for quantitative assessment of severity in pneumoconiosis using both size and frequency of lung nodules that detected by thin-section CT images. This method consists of three steps. First, thoracic organs (body, ribs, spine, trachea, bronchi, lungs, heart, and pulmonary blood vessels) are segmented. Second, lung nodules that have radius over 1.5mm are detected. These steps used functions of our developed computer aided detection system of chest CT images. Third, severity in pneumoconiosis is quantified using size and frequency of lung nodules. This method was applied to nine pneumoconiosis patients. The initial results showed that proposed method can assess severity in pneumoconiosis quantitatively. This paper demonstrates effectiveness of our method in diagnosis and prognosis of pneumoconiosis in CT screening

Visualization and unsupervised clustering of emphysema progression using t-SNE analysis of longitudinal CT images and SNPs

Chronic obstructive pulmonary disease (COPD) is predicted to become the third leading cause of death worldwide by 2030. A longitudinal study using CT scans of COPD is useful to assess the changes in structural abnormalities. In this study, we performed visualization and unsupervised clustering of emphysema progression using t-distributed stochastic neighbor embedding (t-SNE) analysis of longitudinal CT images, smoking history, and SNPs. The procedure of this analysis is as follows: (1) automatic segmentation of lung lobes using 3D U-Net, (2) quantitative image analysis of emphysema progression in lung lobes, and (3) visualization and unsupervised clustering of emphysema progression using t-SNE. Nine explanatory variables were used for the clustering: genotypes at two SNPs (rs13180 and rs3923564), smoking history (smoking years, number of cigarettes per day, pack-year), and LAV distribution (LAV size and density in upper lobes, LAV size, and density in lower lobes). The objective variable was emphysema progression which was defined as the annual change in low attenuation volume (LAV%/year) using linear regression. The nine-dimensional space was transformed to two-dimensional space by t-SNE, and divided into three clusters by Gaussian mixture model. This method was applied to 37 smokers with 68.2 pack-years and 97 past smokers with 51.1 pack-years. The results demonstrated that this method could be effective for quantitative assessment of emphysema progression by SNPs, smoking history, and imaging features

Automated Assessment of Aortic and Main Pulmonary Arterial Diameters using Model-Based Blood Vessel Segmentation for Predicting Chronic Thromboembolic Pulmonary Hypertension in Low-Dose CT Lung Screening

Chronic thromboembolic pulmonary hypertension (CTEPH) is characterized by obstruction of the pulmonary vasculature by residual organized thrombi. A morphological abnormality inside mediastinum of CTEPH patient is enlargement of pulmonary artery. This paper presents an automated assessment of aortic and main pulmonary arterial diameters for predicting CTEPH in low-dose CT lung screening. The distinctive feature of our method is to segment aorta and main pulmonary artery using both of prior probability and vascular direction which were estimated from mediastinal vascular region using principal curvatures of four-dimensional hyper surface. The method was applied to two datasets, 64 low-dose CT scans of lung cancer screening and 19 normal-dose CT scans of CTEPH patients through the training phase with 121 low-dose CT scans. This paper demonstrates effectiveness of our method for predicting CTEPH in low-dose CT screening

Segmentation of aorta and main pulmonary artery of non-contrast CT images using U-Net for chronic thromboembolic pulmonary hypertension : evaluation of robustness to contacts with blood vessels

Enlargement of the pulmonary artery is a morphological abnormality of pulmonary hypertension patients. Diameters of the aorta and main pulmonary artery (MPA) are useful for predicting the presence of pulmonary hypertension. A major problem in the automatic segmentation of the aorta and MPA from non-contrast CT images is the invisible boundary caused by contact with blood vessels. In this study, we applied U-Net to the segmentation of the aorta and MPA from non-contrast CT images for normal and chronic thromboembolic pulmonary hypertension (CTEPH) cases and evaluated the robustness to the contacts between blood vessels. Our approach of the segmentation consists of three steps: (1) detection of trachea branch point, (2) cropping region of interest centered to the trachea branch point, and (3) segmentation of the aorta and MPA using U-Net. The segmentation performances were compared in seven methods: 2D U-Net, 2D U-Net with pre-trained VGG-16 encoder, 2D U-Net with pre-trained VGG-19 encoder, 2D Attention U-Net, 3D U-Net, an ensemble method of them, and our conventional method. The aorta and MPA segmentation methods using these U-Net achieved higher performance than a conventional method. Although the contact boundaries of blood vessels caused lower performance compared with the non-contact boundaries, the mean boundary distances were below about one pixel

Association analysis of SNPs with CT image-based phenotype of emphysema progression in heavy smokers

Chronic obstructive pulmonary disease (COPD) is predicted to become the third leading cause of death worldwide by 2030. Smoking is a well-known risk factor in the development of COPD. Association between COPD genes and smoking have been studied. This paper presents an association analysis of single nucleotide polymorphisms (SNPs) with a CT image-based phenotype of emphysema progression in heavy smokers. The emphysema progression was quantitatively represented by the annual increment of low attenuation volume (LAV) on CT scans for five years. 10 candidate SNPs were selected from 316 SNPs in 125 papers of genetic studies of COPD and lung cancer. The genotypes were determined by real-time polymerase chain reaction (PCR) using deoxyribonucleic acid (DNA) extracted from saliva samples. The association analysis was performed by Fisher's exact test and logistic regression analysis. This method was applied to a dataset with 144 participants (71 smokers, 61 past smokers, and 12 non-smokers). The results showed that the genotypes of rs3923564 and rs13180 SNPs were candidate SNPs associated with the CT image based-emphysema progression

Annual change in bone mineral density in COPD

Background: Osteoporosis is a well-known comorbidity in COPD. It is associated with poor health status and prognosis. Although the exact pathomechanisms are unclear, osteoporosis is suggested to be either a comorbidity due to shared risk factors with COPD or a systematic effect of COPD with a cause–effect relationship. This study aimed to evaluate whether progression of osteoporosis is synchronized with that of COPD. Materials and methods: Data from 103 patients with COPD included in the Hokkaido COPD cohort study were analyzed. Computed tomography (CT) attenuation values of thoracic vertebrae 4, 7, and 10 were measured using custom software, and the average value (average bone density; ABD4,7,10) was calculated. The percentage of low attenuation volume (LAV%) for each patient was also calculated for evaluation of emphysematous lesions. Annual change in thoracic vertebral CT attenuation, which is strongly correlated with dual-energy X-ray absorptiometry-measured bone mineral density, was compared with that in FEV1.0 or emphysematous lesions. Results: In the first CT data set, ABD4,7,10 was significantly correlated with age (ρ=–0.331; p=0.0006), body mass index (BMI; ρ=0.246; p=0.0136), St George’s Respiratory Questionnaire (SGRQ) activity score (ρ=–0.248; p=0.0115), eosinophil count (ρ=0.229; p=0.0198), and LAV% (ρ=–0.372; p=0.0001). However, ABD4,7,10 was not associated with FEV1.0. After adjustment for age, BMI, SGRQ activity score, and eosinophil count, no significant relationship was found between ABD4,7,10 and LAV%. Annual change in ABD4,7,10 was not associated with annual change in LAV% or FEV1.0. Conclusion: Progression of osteoporosis and that of COPD are not directly related or synchronized with each other

Automated detection method of thoracic aorta calcification from non-contrast CT images using mediastinal anatomical label map

Progression of thoracic aortic calcification (TAC) has been shown to be associated with hard cardiovascular events including stroke and all-cause mortality as well as coronary events. In this study, we propose an automated detection method of TACs of non-contrast CT images using mediastinal anatomical label map. This method consists of two steps: (1) the construction of a mediastinal anatomical label map, and (2) the detection of TACs using the intensity and the mediastinal anatomical label map. The proposed method was applied to two non-contrast CT image datasets: 24 cases of chronic thromboembolic pulmonary hypertension (CTEPH) and 100 non-CTEPH cases of low-dose CT screening. The method was compared with two-dimensional U-Nets and the Swin UNETR. The results showed that the method achieved significantly higher F1 score of 0.937 than other methods for the non-CTEPH case dataset (p-value < 0.05, pairwise Wilcoxon signed rank test with Bonferroni correction)

An automated distinction of DICOM image for lung cancer CAD system

Automated distinction of medical images is an important preprocessing in Computer-Aided Diagnosis (CAD) systems. The CAD systems have been developed using medical image sets with specific scan conditions and body parts. However, varied examinations are performed in medical sites. The specification of the examination is contained into DICOM textual meta information. Most DICOM textual meta information can be considered reliable, however the body part information cannot always be considered reliable. In this paper, we describe an automated distinction of DICOM images as a preprocessing for lung cancer CAD system. Our approach uses DICOM textual meta information and low cost image processing. Firstly, the textual meta information such as scan conditions of DICOM image is distinguished. Secondly, the DICOM image is set to distinguish the body parts which are identified by image processing. The identification of body parts is based on anatomical structure which is represented by features of three regions, body tissue, bone, and air. The method is effective to the practical use of lung cancer CAD system in medical sites